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ABSTRACT 



This is the first paper of a series in which we will attempt to put constraints on the flattening of dark halos in disk galaxies. We observe 
for this purpose the HI in edge-on galaxies, where it is in principle possible to measure the force field in the halo vertically and radially 
from gas layer flaring and rotation curve decomposition respectively. In this paper, we define a sample of 8 HI rich late-type galaxies 
suitable for this purpose and present the HI observations. 
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1. Introduction 

Since the 1970's, it has been known that the curvature of the 
universe is remarkably flat. This implies that the ratio of the to- 
tal density to critical density is Q. to , « 1. At this time, it was 
also known from measurements of stars and gas that the mass 
density of luminous matter is Q/,„„ < 0.005 — less than 0.5% of 
that required for a flat Sl m ~ 1 universe. Indeed , the missing mat- 
ter con troversy had be gun in the 1930's w hen lOortl ( 1 19321) . fol- 
lowing lKaptevnl (1 1 9221) . and lZwickvld 19371) independently found 
evidence for vast amounts of unseen matter on different scales. 
Oort's analysis of the spatial and velocity distribution of stars 
in the Solar neighbourhood concluded that luminous stars com- 
prised approximately half the total mass indicated by their mo- 
tion, assuming gravititational equilibrium. Zwicky's analysis of 
the velocity dispersions of rich clusters found that approximately 
90 - 99% of the mass was unseen, if the systems were gravita- 
tionally bound. Later on, the inventory of the mass in the Solar 
neighbourhood showed that the dark matter in the galactic disk 
at the solar galact ocentric radius was mostly fainter stars and in- 
terstellar gas (e.g. lHolmberg & FlvnnL 120001) . Zwicky's value of 
the velocity dispersion in the Coma cluster is close to the cur- 
rent one (e.g. lColless & Dunnlll996l) . and, despite the presence 
of hot X-ray emitting gas, there must be dark matter to cause 

^matter ~ 0.2 — 0.3. 

Modern dark matter research began in 1970 with several pa- 
pers which found that galaxies contain more grav i tating matter 
than can be accounted for by the stars. iFreemanl (1 19701) noted 
that the atomic hydrogen (HI) rotation curves of the late-type 
disk galaxies NGC300 and M33 peaked at a larger radius than 
expected from the stellar light distribution. This implied that 
the dark matter was more extended than the stellar distribu- 



tion and that its mass was at least as great as the luminous 
mass. Howeve r, these data were of poor angular resolution. At 
the same time iRubin & Fordl (1 19701) published a rotation curve 
of M31 based on optical data, which did not seem to decline 
in the outer parts. In the 1970s radio observations of increas- 
ingly better angular resolution and better sensitivity showed 
that flat rotation curves are typical in spiral galaxies ( Shostak , 
197l HRogstad & ShostakLll972tlRoberts & RotsL [l973t iBosma , 



19781) and that the HI extent of a spiral galaxy can be far greater 
than the extent of the optical image. This, combined with surface 
photometry, leads to very high mass-to-light ratios in the o uter 
parts, as shown clearly for seve ral galaxies in lBosmal d 19781) and 
Bosma & van der Kruitl (119791) . 

In 197 4, dark matter halos were found to extend eve n fur- 
ther, when lOstriker et al.l d 19741) and lEinasto et al.1 d 19741) tabu- 
lated galaxy masses as a function of radius and found that galaxy 
mass increased linearly out to at least 100 kpc and 10 12 M for 
normal spirals and ellipticals. Despite this large dark matter frac- 
tion inferred in galaxies, it was still not suffici ent to reach the 
critica l D., n , ~ 1 value of a flat universe. See e.g. ISofue & Rubinl 
(1200 ll) and lRobertsI d2008h for a brief history of dark matter in 
galaxies. 

Around the same time, application of nuclear physics to Big 
Bang theory showed that big bang nucleosynthesis (BBNS) in 
the early universe produced specific abundances of the light ele- 
ments, and predicted the total baryon density was =s 0.044. In 
the last decade, the combined observations of high redshift SNIa, 
dReiss et al. L 119981: iPerlmutter etaiT. 1 19991). the 2-Degree Field 
Galaxy Redshift Survey (2dFGRS) dPercival et all 1200 lb and 
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1 KCF and PCvdK recall influential colloquia and other presentations 
by M.S. Roberts in the early 1970s on a flat HI rotation curve of the 
Andromeda galaxy which helped to steer the evolution of this subject. 
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Table 1. HI observations 
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WMA P microwave background measurements (ISpergel et all 
120031) confirmed the baryonic mass density O^, found by BBNS 
and established that dark energy comprises about 75% of the 
critical density. Consequently, the scale of the missing dark mat- 
ter is now reduced. The mass density Q m is now only 0.25, but 
the problem remains: the measured baryonic mass density is still 
only ~18% of the total mass density Q m . Thus, dark matter ac- 
counts for 82% of mass in the universe, and baryonic matter is 
only 4.5% of the total content of the universe^ 

More than 1000 galactic rotation curves have now been mea- 
sured, and very few display a Keplerian decline with radius. HI 
and Ha rotation curves of spiral galaxies show that the total- 
mass-to-light ratios are typically M/L = 10-20 M /L , and the 
luminous matter therefore accounts for only 5 - 10% of the 
total mass inferred from the rot ation curves. F or low surface 
brightness (LSB) dis k galaxies (Ide Blokl 1 19971) and dwarf ir- 
regular (dl) galaxies (ISwatersL 1 19991) . the M/L values increase 
to 10-100 M /L , with an extreme of 220 for ESO215-G009, 
a gas-rich LSB galaxy with a very high gas mas s to light ratio 
of M ini Lb = 21 and low recent star formation dWarren et al.L 
12004 . 

Dwarf spheroidal (dSph) galaxies, with typical total masses 
of only ~10 7 M within radii of a few hundred parsecs, are even 



2 For this illustrati on we used the pa rameters adopted in the 
Millenium Simulation (Springel et al., 2005), which are fl,„ = Q. dm + 
Sl b = 0.25, Sl b = 0.045, n A = 0.75. 



more extreme: several have very large dark matter fractions with 
mass-to-light ratios in the range 200 - 1000 M /L . In these 
faint, small galaxies, the dynamical mass is estimated from the 
line-of-sight ve l ocities of individual stars. The Ursa Major dSph 
dWillman et all 120051) , is one of the most dark matter domi- 
nated galaxies known to date with a central mass-to-light ratio 
M/L ~500 M /L , which is believed to increase further at larger 
radii dKlevna e t al., 2005). These systems appear to have only 
very small baryonic mass fractions. 

The rotation curves of disk galaxies are important probes of 
the equatorial halo potential gradient. By decomposing the ob- 
served rotation curve into contributions from the visible mass 
components, the radial potential gradient of the halo can be 
measured, assuming the system is in centrifugal equilibrium. 
The dark halo mass component is typically fitted by a pseudo- 
isothermal halo model with density 



p(R)=Po 



1+1* 

Rr 



(1) 



where the halo is characterised by its central density po and core 
radius R c . Pseudo-isothermal halos have an asymptotic density 
p oc R 2 at large radii which is consistent with commonly flat 

rot ation curves. 

iKormendv & Freemanl (l2004f) compiled the published dark 
halo density distributions for a large sample of Sc-Im and dSph 
galaxies, and found well-defined scaling relationships for the 
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Table 2. Resolution of HI observations 



Table 3. Noise of HI channel maps 
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dark halo parameters for galaxies spanning a range of 6 decades 
of luminosityQ They found that halos of less luminous and 
massive dwarf spheroidals have higher central densities, up to 
~1 M pc -3 and core radii of ~0.1 kpc, compared to ~ 10~ 3 
M pc~ 3 and ~30-100 kpc respectively, for large bright Sc 
galaxies. The observed correlations suggest a continuous phys- 
ical sequence of dark halo population in which the properties 
of the underlying dark halo scale with the baryon luminosity 
(iKormendv & FreemanLl2004 . 

We now consider the flattening of the dark halo density dis- 
tribution, defined by q — c/a, where c is the halo polar axis and 
a is the major axis in the galactic plane. The vertical distribution 
of the halo is much more difficult to measure than the radial dis- 
tribution in the equatorial plane, as most luminous tracers of the 
galaxy potential gradient lie in the plane of the galaxy and offer 
little indication of the vertical gradient of the potential. Past mea- 
surements that were obtained with a variety of different methods 
gave a large range of q from 0.1 to 1.37, with no concentration 
on any particular value. For our Galaxy, the halo shape has been 
measured more than ten times by four different methods, that 
yielded ^-values ranging from 0.45 to 1.37. 

In this series of papers, we will use the flaring of the HI gas 
layer to measure the vertical flattening of the dark halo, because 
we believe this method to be the most promising for late-type 
disk galaxies. Like all steady state mass components of a galaxy, 
the gaseous disk of an isolated disk galaxy can be assumed to 
be in hydrostatic equilibrium in the gravitational potential of the 
galaxy, unless there are visible signs from the gas distribution 
and kinematics that the HI layer is disturbed (e.g. by mergers 
or local starbursts). In the vertical direction, the gradient of gas 
pressure is balanced by the gravitation force (ignoring any con- 
tribution from a magnetic pressure gradient). From the observed 
distribution of the gas velocity dispersion and the gas density 
distribution, we can in principle measure the total vertical gravi- 
tational force. 

Euler's equation for a steady-state fluid of density p, veloc- 
ity v and pressure p in a gravitational potential <D is -(v.V)v = 
p _1 Vp + V<D. In the case of a vertically Gaussian gas density dis- 
tribution with a vertically isothermal gas velocity dispersion, the 
gradient of the total vertical force K z in the z-direction can be 
calculated directly from the gas velocity dispersion cr v g (R) and 
the gas layer thickness FWH1VL g (7?), each of which are functions 
of radius: 



dKz 
dz 



(FWHM ; , g /2.35) 2 



(2) 



3 E-Sbc galaxies were not included to avoid the larger uncertain- 
ties associated with stellar bulge-disk decomposition and the relatively 
larger contribution of the stellar mass of varying stellar ages. 



It is possible to measure the halo shape over the entire HI ex- 
tent of the luminous disk using the flaring of the HI distribution, 
which typically extends in radius to 2 - 3R 25, and by measuring 
the density distribution of the gas and stellar distributions. For 
a given vertical gas velocity dispersion, a more flattened dark 
halo requires decreased flaring and increased gas surface den- 
sity. The relatively high gas content of late-type galaxies allows 
measurement of both the halo vertical force field from the gas 
layer flaring and the halo radial force field from rotation curve 
decomposition. In this series of papers we will attempt to mea- 
sure the halo flattening using the flaring of HI disk in eight small 
late-type disk galaxies. 

The main advantage of this method is that it can be used for 
all gas-rich spiral galaxies inclined sufficiently to measure accu- 
rate kinematics, unlike some methods which are applicable only 
to specific kinds of galaxies like polar ring galax ies. Th i s mini - 
mum inclination was determi ned to be ~ 60° bylOllind d 19951) . 
This method was first tried bv lCelnik et al.l(l 19791) on the Galaxy , 
and early development was undertaken by Ivan der Kruitl (Il98ll) . 
who applied it to low resolution observations of NGC891 and 
concluded that the halo was not as flattened as the stellar disk. 
It was then used for several galaxies in the 1990's, most notably 
the careful study of the very nearby Sc galaxy NGC4244 which 
was found to ha ve a highly fla ttened halo with q = 0.2^ 3 out to 
radii of * 2R 25 (lOllin Jll996h . 

All applications of the HI flaring method to date have de- 
rived highly flattened halo distributions (q < 0.5). With the 
exclusion of NGC4244, we suspect that the assumption of a 
radially-constant gas velocity dispersion led to errors in the mea- 
sured total vertical force, and thereby to the derived flatten- 
ing of the halo. The other difficulty for all these galaxies ex- 
cept NGC4244 is that they are large galaxies with peak rotation 
speeds > 170 km s _1 . Given that ( i) spiral galaxies typically have 
HI velocity dispersions in the relatively small range 6-10 km 
s~ 1 (lTamburro et all 120091) and (ii) the halo shape q is roughly 
proportional to the gradient of the vertical force, we see from 
Eqn. |2] that the HI thickness is to first order inversely propor- 
tional to the peak rotation speed. Consequently, the HI flaring 
method should work best for small disk galaxies with relatively 
low rot ation speeds, as these galaxies should exhibit more HI 
flaring. iBosmal (1 19941) already showed by calculations that flar- 
ing is relatively more important in galaxies with low circular 
velocity. 

The flaring method ha s also been applied to the Galaxy by 
lOlhng & Merrifieidl(l2000l) . However, uncertainty in the position 
and rotation speed of the Sun resulted in a large uncertainty of 
the measured halo shape 0.5 < q < 1.25. 

In this first paper we discuss the selection of our sample and 
present HI observations. In paper II we will review methods to 
derive the information we need for our analysis from these data, 
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Table 4. HI measurements 
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namely the rotation curves and the HI distributions, HI velocity 
dispersion and the flaring of the HI layer, all as a function of 
radius in the deprojected galaxy plane. Paper III will be devoted 
to applying this to our data and presenting the results for each 
individual galaxy. In paper IV will analyse the data of one galaxy 
in our sample, namely UCG7321 and we will set limits on the 
flattening of its dark halo. 

2. Observations 

2. 1 . Sample selection 

We selected a sample of southern disk galaxies, for observation 
with the Australia Telescope Compact Array (ATCA), that were 
small, Hl-rich and late-type, in order to maximise the expected 
flaring and the likelihood that the HI brightness would be suffi- 
cient to probe the low surface densities needed to measure the 
vertical structure of disk galaxies. 

The galaxies were chosen to be edge-on (optical major-to- 
minor axis ratio a/b > 10) to simplify HI modelling and increase 
the accuracy of the measured HI density and kinematics. We 
avoided galaxies at Galactic latitude \b\ < 10° to reduce Galactic 
extinction in the optical bands. To minimise uncertainties in the 
stellar mass distribution we chose relatively bulgeless galaxies 
with Hubble type Scd-Sd. The galaxies were chosen to be nearer 
than 30 Mpc to enable the flaring of the HI to be resolved with 
the 6" ATCA synthesised beam. At 30 Mpc the FWHM of the 
ATCA beam is ~1 kpc. The vicinity of each galaxy was searched 
in NED for nearby neighbours to ensure that it was isolated. The 
HI masses of all galaxies satisfying these criteria were measured 
from the HI Parkes All-Sky survey (HIPASS) online data re- 
leastfl With a minimum HI flux integral of 15 Jy km s _1 , this 
limited our sample to only 5 galaxies. 

The nearest isolated thin edge-on galaxy, ESO274-G001, 
was added to the sample. Although it had a Galactic latitude 
of 9.3°, it is exceptionally close at 3.4 Mpc in the Centaurus 
A galaxy group, allowing the HI disk to be measured at a resolu- 
tion of ~ 150 pc. A search of HI observations in the VLA archive 
revealed the superthin northern Sd galaxy UGC7321 which we 
also added to our sample. We also added the Sb galaxy IC2531 



-69" 18' — 




Right Ascension 



Fig. la. ESO074-G015 (IC5052). HI column density map 
overlaid on the DSS image. Contours are plotted in 
(3,5, 10,25,50,80) x a, where the rms noise cr = 2.05 mJy 
beairr'or 2.79x 10 19 atoms cirT 2 . The FWHM synthesised beam 
has dimensions 9.0" x 9.0", and is displayed as a cross-hatched 
symbol in the lower left. 



http://www.atnf.csiro.au/research/multibeam/release 



to our sample, since it has an extended HI disk and a large quan- 
tity of archival observations available. The range of maximum 
rotation velocity of our sample galaxies in shown in column 
7 of Table |4] With the inclusion o f IC2531, a galaxy with a 
box/peanut bulge, and hence barred dBureau &FreemanLll999l) . 
we attempt to measure the halo shape of galaxies spanning dif- 
ferent mass scales and stages of secular evolution. 
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Fig. lb. ESO074-G015 (IC5052). Top left: XV map. XV con- 
tours are (3,5, 10,25,50, 100)x<x, where therms noise cr = 1.90 
mJy bearrr'or 2.58 x 10 19 atoms cirr 2 . The half power beam 
extent over the major axis is shown in the lower left corner. 
Top right: Integrated spectrum. Middle left: HI column den- 
sity map rotated with the galaxy major axis aligned with the X 
axis. Column density map contours are (3, 5, 10, 25, 50, 100)x<x, 
where the rms noise cr - 1.65 mJy beam -1 or 2.24 x 10 19 atoms 
cm -2 . The synthesised beam has dimensions 9.0" x 9.0", and 
is displayed as a cross-hatched symbol in the lower left. Middle 
right: Minor axis profile. Bottom left: Major axis profile. 



2.2. Observing modes and method 

Galaxies were observed with the ATCA in a range of array con- 
figurations to obtain high spatial resolution across each galaxy. 
In addition, available observations from the ATCA and VLA 
archives were also used. The list of observations is shown in 
Table Q] Column 8 shows the integration time in hours of each 
observation. The resulting resolution and sensitivity of the HI 
observations are shown in Table [2] and Table [3] respectively. 
We have smoothed the data in Right Ascension to match that in 
Declination, providing the smallest circular beam possible. The 
spatial resolutions shown in Table [2] correspond to the FWHM 
of the resulting beams. 

At the ATCA the XX and YY polarisations were used, with 
a spectral channel width of 3.3 km s . At the VLA the RR and 
LL polarizations were used, and the channel width was 5.2 km 
s _1 . Due to the narrow bandwidth available with the VLA cor- 
relator at the time of observations, multiple overlapping band- 
passes were observed to fully span the velocity width of the 
target galaxies. For those observations, each bandpass was cali- 
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Fig. lc. ESO074-G015 (IC5052). Top : Peak brightness temper- 
ature map. Contours are plotted in (35,75, 120) x cr, where the 
rms noise is 10.0 K. The FWHM synthesised beam has dimen- 
sions 9.0" x 9.0". Middle: Major axis peak brightness temper- 
ature profile. Bottom: Inferred HI opacity calculated assuming 
constant HI spin temperatures. The resulting maximum opaci- 
ties along each line of sight column through the galaxy disk are 
plotted for T spin = 200, 300 and 400 K (bottom to top). 



brated and imaged separately, and the resulting subcubes glued 
together along the spectral axis. 



3. Data reduction and imaging 

A detailed account and discussion of the data reduction and 
imaging procedures is given in the online appendix. Here we 
restrict ourselves to a brief synopsis. 

The reduction of the HI data was performed with the radio 
interferometry data reduction package miriad. The data were cal- 
ibrated in the usual manner, using primary and secondary flux 
calibrators. Solar and terrestrial interference peaks in the cali- 
brator scans were inspected, and flagged if necessary. The band- 
pass calibration was done using either the primary or, in case 
of sufficient signal-to-noise ratio, the secondary calibrator. After 
calibration, the target data were inspected for interference, and 
flagged appropriately. The radio continuum contribution to each 
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Fig. 2a. ESO109-G021 (IC5249). HI column density map over- 
laid on the DSS image. Contours are plotted in (3, 5, 10, 25, 50)x 
cr, where the rms noise cr = 1.81 mJybeam _1 or3.11xl0 19 atoms 
cm -2 . The FWHM synthesised beam has dimensions 8 .0" x 8 .0", 
and is displayed as a cross-hatched symbol in the lower left. 



channel map was removed in the uv-plane, by using a suitable 
average of line-free channel maps. 

4. Hi emission 

4. 1 . The HI distribution 

To make the HI column density map, the intensity of each chan- 
nel map was first converted from flux density I(x, y) in units of 
Jy beam -1 to HI column density Nm in atoms crrT 2 using 

M /1.101974351 xl0 24 \ 
tffflOO = I(x,y), 



(3) 



where 6 ma j and #„„„ are the major and minor axis FWHM of the 
synthesised beam in arcsec. The HI column density maps were 
made by taking the zeroth moment of masked channel maps 
(with v in km s _1 ): 



■I 



N m (v)dv. 



(4) 



To exclude any residual sidelobes, a loose region mask was 
defined interactively for each channel map and then all pixels 
with flux density above a nominal clip level of 2-3<x in each 
masked channel map were integrated. This procedure ignores 
faint, low-level HI emission and the resulting integrated HI maps 
will not contain any extended, low surface brightness component 
in the gas distribution. However, this does not have an effect on 
the general distribution of HI as determined by our method. The 
maps in this paper have been produced to show the general dis- 
tribution of the gas in the disk in order to judge the suitability 
of the galaxies for further study. We will only be able to anal- 
yse the flaring and velocity dispersion as long as we can obtain 
high signal-to-noise profiles. A search for extended HI is beyond 
the scope of our present analysis of the data. As the number of 
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Fig. 2b. ESO109-G021 (IC5249). Top left: XV map. XV con- 
tours are (3, 5, 10, 20, 40) x cr, where the rms noise cr - 1.76 mJy 
bearrr'or 3.03 x 10 19 atoms crrT 2 . The half power beam extent 
over the major axis is shown in the lower left corner. Top right: 
Integrated spectrum. Middle left: HI column density map rotated 
with the galaxy major axis aligned with the X axis. Column den- 
sity map contours are (3, 5, 10, 25, 50) x cr, where the rms noise 
o- = 1.62 mJy beam" 'or 2.79 x 10 19 atoms cm" 2 . The FWHM 
synthesised beam has dimensions 8.0" x 8.0", and is displayed 
as a cross-hatched symbol in the lower left. Middle right: Minor 
axis profile. Bottom left: Major axis profile. 



channels integrated to form the column density map varies as 
a function of position, the noise in the HI column density map 
also varies accordingly. In order to approximate the noise of the 
HI column density map, a noise map was also constructed. The 
noise at each position in the map was calculated as 



CTtotal = V"^h X <CT ch (v)>, 



(5) 



where (cr ch (v)) is the mean rms noise of all channel maps inte- 
grated at that position. By dividing the column density map by 
the noise map a signal-to-noise rat io map was formed. Following 
the method of iKregel etal] d2004l) . the rms noise of the column 
density map was calculated from the average flux density of all 
pixels with a signal-to-noise ratio between 2.75 and 3.25. 



4.2. Global properties 

The integrated HI spectrum was measured from the HI data cube 
using the HI column density map to mask the cube, and integrat- 
ing over the area of the galaxy mask. This flux density spectrum 
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Fig.2c. ESO109-G021 (IC5249). Top : Peak brightness temper- 
ature map. Contours are plotted in (35,58,93) x <x, where the 
rms noise is 11.7 K. The FWHM synthesised beam has dimen- 
sions 8.0" x 8.0". Middle: Major axis peak brightness temper- 
ature profile. Bottom: Inferred HI opacity calculated assuming 
constant HI spin temperatures. The resulting maximum opaci- 
ties along each line of sight column through the galaxy disk are 
plotted for T spi „ = 200, 300 and 400 K (bottom to top). 



was converted to units of Jy by dividing by the synthesised beam 
area 



S(v) = 



1 



1.133 6 ma j @min 



y)dxdy. 



(6) 



As the flux was integrated over an equal area in each channel 
map, the noise of each channel is 



CTto 



-\/Area 



hi x <r ch , 



(7) 



where Areani is the area of HI column density emission, and <x ch 
is the rms noise of each channel map. For HI cubes with a con- 
stant rms noise in each channel, the noise of the spectrum was 
fitted from the line-free channels in the spectrum. The velocity 
widths at the 20%- and 50%-levels (W20, W50, respectively) were 
measured from the spectrum, and the spectrum was integrated 



Fig. 3a. ESO115-G021. HI column density map overlaid on the 
DSS image. Contours are plotted in (3, 5, 10, 25, 50) x <x, where 
the rms noise cr — 1.76 mJy beairr'or 2.44 x 10 19 atoms cm' 2 . 
The FWHM synthesised beam has dimensions 8.9" x 8.9", and 
is displayed as a cross-hatched symbol in the lower left. 



over the velocity range spanned by W20 to get the integrated flux, 
FT: 



Fl 



JW20.I, 



W20Mgh 



S(v)dv. 



(8) 



The restriction of the integration to the velocity range W20 in- 
troduces a small but systematic error in the total flux, underes- 
timating it by no more than a few %. However, for our present 
purposes it is important to examine the integrated profile in or- 
der to identify major asymmetries that would make the galaxy 
unsuitable for our purposes. 

The total HI mass, Mui was then measured from the flux 
integral using the formula 



M m = 2.35D 2 x FI, 



(9) 



where the adopted distance D is given in Mpc. These HI mea- 
surements are shown in Table |4] These total masses are also un- 
derestimated by a few %. 

The centre of each galaxy was obtained from the HI column 
density map by rotating the image by 180° and finding the opti- 
mal pixel offset of the rotated image relative to the original image 
by maximising the correlation function of the two images using 
the IDL Astronomy User's Library function correLoptimizeu 
Formally this gives the center of the HI disk rather than the 
galactic rotational center. Comparison to the coordinates given 
in the NASA/IPAC Extragalactic Database NEE0, which are de- 
rived from the Two Micron All Sky Survey 2MASS0, show dif- 
ferences of on average only 4 arcsec per coordinate, except for 
ESO074-G015, where it is about 20 arcsec in RA and 30 arc- 
sec in Dec. However, this galaxy shows some deviations from 



http://idlastro.gsfc. nasa.gov/homepage.html. 

6 nedwww.ipac.caltech.edu. 

7 
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Fig.3b. ESO115-G021. Top left: XV map. XV contours are 
(3,5, 10,20,50, 100, 150)xcr, where therms noise cr = 1.55 mJy 
bearrr'or 2.16 x 10 19 atoms cm -2 . The half power beam extent 
over the major axis is shown in the lower left corner. Top right: 
Integrated spectrum. Middle left: HI column density map rotated 
with the galaxy major axis aligned with the X axis. Column den- 
sity map contours are (3, 5, 10, 25, 50) x cr, where the rms noise 
cr = 1.47 mJy beam" 1 or 2.04 x 10 19 atoms cnr 2 . The FWHM 
synthesised beam has dimensions 8.9" x 8.9", and is displayed 
as a cross-hatched symbol in the lower left. Middle right: Minor 
axis profile. Bottom left: Major axis profile. 



symmetry, both in optical appearance as well as in the HI distri- 
bution. The position adopted fits the symmetry of the XV map 
very well (see Fig. [TBI and is also closer to centers listed in NED 
that are derived from photographic and IRAS data. 

The position angle of the galaxy was found by a similar 
method. First a wide range of position angles were trialed around 
the estimated position angle. The HI column density map was re- 
flected over each test position angle about the galaxy center and 
the residual of the two maps was calculated. The position an- 
gle that yielded the lowest summed residual was then adopted 
as the next position angle estimate. Subsequent rounds spanned 
smaller and smaller ranges of the position angle until the posi- 
tion angle was determined to less than 0.1 deg or until the to- 
tal residual dropped below 1% of the summed flux of the HI 
column density map. Both methods worked very well for ob- 
servations with a near circular synthesised beam, achieving an 
accuracy of ~1 pixel in the centre position and 0.1 - 0.2 deg in 
position angle. However, observations with a highly elongated 
beam were poorly fitted, as the distortion of the image caused 
by the beam shape biased the derived PA towards the PA of the 




Fig. 3c. ESO115-G021. Top : Peak brightness temperature map. 
Contours are plotted in (35,42,86,129) x cr, where the rms 
noise is 8.6 K. The FWHM synthesised beam has dimensions 
8.9"x8.9". Middle: Major axis peak brightness temperature pro- 
file. Bottom: Inferred HI opacity calculated assuming constant 
HI spin temperatures. The resulting maximum opacities along 
each line of sight column through the galaxy disk are plotted for 
T spi „ = 200, 300 and 400 K (bottom to top). 



synthesised beam. For such observations the uncertainty of the 
measured centre and position angle were larger. 

The xzv cube was constructed by rotating the HI data cube 
about the galaxy center to align the galaxy major axis with the 
X axis. In Sect. [5] we show the rotated channel maps to allow 
the flaring to be viewed directly from the observations. We also 
present the HI column density map with the projected major 
and minor axis profiles as measured from the HI column den- 
sity map. The HI diameter and maximum vertical extent were 
measured at the 3<x level of the projected profiles. 

An XV cube was also produced by reordering the cube axes 
to form a xvz cube. This cube was then integrated over the z 
axis to make an XV map using the same method as was used to 
make the HI column density map. Both the XV map and the HI 
column density map are shown with their respective noise levels 
in Sect. [5] The systemic velocity of the galaxy was measured 
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Fig. 4a. ESO138-G014. HI column density map overlaid on the 
DSS image. Contours are plotted in (3,5, 10,25,50, 100) x a, 
where the rms noise <x = 2.84 mJy bearrr'or 2.73 x 10 19 atoms 
cm -2 . The FWHM synthesised beam has dimensions 10.7" x 
10.7", and is displayed as a cross-hatched symbol in the lower 
left. 



from the XV map by shifting the XV distribution flipped in V 
with respect to the actual XV map to maximise the correlation 
function. 

Peak brightness temperature maps were also made for each 
galaxy, where the peak brightness temperature is the maximum 
brightness temperature over all channel maps at each spatial po- 
sition. The brightness temperature Tg was derived from the HI 
column density Nui in each channel map using 



T R = 



1.83 x 10 18 ' 



(10) 



The major axis profile of the peak brightness map is also shown 
for each galaxy in Sect. [5] These plots are useful indicators of 
possible HI self-absorption that could be obscuring the intrin- 
sic HI surface distribution. Sect. [5] also shows the inferred HI 
opacity of the brightest HI emission at each major axis posi- 
tion, assuming three different values of the HI spin temperature 
(200, 300, 400 K). 



5. Results for individual galaxies 

In this section we present our data for the galaxies in the sample. 
For each galaxy we present the following figures. Part a shows 
an HI column density map overlaid on an optical image from 
the Digitized Sky Survey (DSSfl Part b shows the total inten- 
sity position-velocity map (XV map) together with the integrated 
spectrum, and the HI column density map aligned with the major 



8 The Digitized Sky Surveys were produced at the Space Telescope 
Science Institute under U.S. Government grant NAG W-2166. The im- 
ages of these surveys are based on photographic data obtained using the 
Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt 
Telescope. The plates were processed into the present compressed dig- 
ital form with the permission of these institutions. 
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Fig. 4b. ESO138-G014. Top left: XV map. XV contours are 
(3, 5, 10, 20, 50) x <x, where the rms noise a = 2.74 mJy beam" 1 
or 2.64 x 10 19 atoms cm" 2 . The half power beam extent over the 
major axis is shown in the lower left corner. Top right: Integrated 
spectrum. Middle left: HI column density map rotated with the 
galaxy major axis aligned with the X axis. Column density map 
contours are (3, 5, 10, 25, 50, 75, 100) x cr, where the rms noise 
a = 2.38 mJy beam" 1 or 2.29 x 10 19 atoms cm" 2 . The FWHM 
synthesised beam has dimensions 10.7"x 10.7", and is displayed 
as a cross-hatched symbol in the lower left. Middle right: Minor 
axis profile. Bottom left: Major axis profile. 



axis of the galaxy and the corresponding major and minor axis 
surface density profiles. Part c shows the peak brightness tem- 
perature map and the corresponding observed major axis profile 
and compares this with opacity profiles calculated with some as- 
sumed HI sp in temperatures . The individual channel maps were 
presented bv lO'Brienl d2007l) and are included in further parts d, 
e etc. (depending on how many figures are needed) as online- 
only material at the end of this paper. 



5.1. ESO074-G015 (IC5052) 

ESO074-G015 is a SBd galaxy with a substantial star formation 
region extending above and below the galactic plane on the west 
side of the galactic centre. However despite this, the HI column 
density distribution (shown in Fig. [Tbb is surprisingly symmet- 
ric. Although there is extensive high latitude HI gas extending 
to 3.7 kpc above and 4.2 kpc below the galactic equator, the dis- 
tribution of elevated gas appears to be similar on both sides of 
the galactic centre. Similarly the XV diagram is also fairly sym- 
metric. There is a slight decrease in flux along the line-of-nodes 
of the west side of the galaxy possibly indicating gas depletion 
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Fig. 4c. ESO138-G014. Top : Peak brightness temperture map. 
Countours are plotted in (41,51,83, 124) x cr, where the rms 
noise is 10.4 K. The FWHM synthesised beam has dimensions 
10.7" x 10.7"r. Middle: Major axis peak brightness tempera- 
ture profile. Bottom: Inferred HI opacity calculated assuming 
constant HI spin temperatures. The resulting maximum opaci- 
ties along each line of sight column through the galaxy disk are 
plotted for T spi „ = 200, 300 and 400 K (bottom to top). 



due to a star burst. However there is no apparent localised region 
of high velocity dispersion gas that could be associated with the 
star burst. 

The HI distribution extends to a maximum radial extent of 
281" or 9.1 kpc (for the adopted distance of 6.7 kpc). The scale- 
length of the projected distribution hx is ~ 120 " or 3.9 kpc, with 
the edge of the HI disk at 2.3/ix- The exponential scaleheight of 
the projected minor axis distribution h- is 20.5 " or 667 pc. The 
high latitude HI clouds extend to a height of 3.9 kpc or 6 scale- 
heights away from the equatorial plane. 

The galaxy is well-resolved by the synthesised beam along 
the major axis, extending to 31 FWHM beamwidths on each 
side of the galactic centre, where each FWHM beamwidth is 
292 pc. In the vertical direction the measured scaleheight of the 
projected distribution is approximately 4x HWHM beamwidth, 
indicating that the vertical HI structure is well-resolved by the 
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Fig. 5a. ESO146-G014. HI column density map overlaid on the 
DSS image. Contours are plotted in (3, 5, 10, 25, 50) x cr, where 
the rms noise cr — 1.53 mJy beam -1 or 2.92 x 10 19 atoms cm" 2 . 
The FWHM synthesised beam has dimensions 7.6" x 7.6", and 
is displayed as a cross-hatched symbol in the lower left. 



synthesised beam. Indeed the high latitude HI gas extends to 13x 
HWHM beamwidth. Vertical flaring of the gas distribution can 
be noted in the extreme velocity channels 499, 509, 657, 667, 677 
km s in the outer part of the galaxy where the rotation curve 
inferred from the XV map is roughly flat. 

The HI column density distribution is dominated by two 
bright spots at a projected radius of ~1 kpc from the centre. The 
brightest gas emission in the cube with brightness temperatures 
of approximately 140 K also correspond to this position. The 
HI disk also contains two outer bright points (although less de- 
fined on the western side) at a radius of ~5-6 kpc with a equa- 
torial brightness distribution similar to that of the barred galaxy 
UGC7321, which also displays the characteristic "figure-8" pat- 
tern of a bar in its XV distribution. Although ESO074-G015 is 
classified as a barred SBd galaxy in NED, it does not display a 
"figure-8" pattern in the XV map. ESO074-G015 contains only 
1 /3 of the HI mass of UGC732 1 , and has a lower maximum rota- 
tion speed of 90. 1 km s _1 , relative to 11 2.1 km s" 1 forUGC7321. 
It is possible that the lower angular momentum and overall mass 
of ESO074-G015 restricts gas flow in the orbital patterns of a 
barred potential. 

5.2. ESO109-G021 (IC5249) 

ESO109-G021 is a superthin Sd galaxy with B-band stel- 
lar sc alelength-to-scaleheight ratio of 11 (Ivan der Kruit et ail 
l200~hV The rota tion curve was i nitially fo und to be in solid 
body r otation by lAbe et al.l (1 19991) . however Ivan der Kruit et al] 
(1200 lh pointed out such a rotation curve was inconsistent with 
the double-h orned HI spectrum ( see Fi g.l2bl. By modelling the 
XV diagram lyander Kruit et all j200Tb show that ESO109-G21 
is differentially rotating. 

To improve signal-to-noise and recovery of image structure 
we obtained two additional 12-hour synthesis observations at the 
ATCA and added these to the earlier observation by Carignan 
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Fig. 5b. ESO146-G014. Top left: XV map. XV contours are 
(3,5, 10, 25, 50) x cr, where the rms noise cr — 1 .45 mJy beam -1 
or 2.76 x 10 19 atoms cm 2 . The half power beam extent over the 
major axis is shown in the lower left corner. Top right: Integrated 
spectrum. Middle left: HI column density map rotated with the 
galaxy major axis aligned with the X axis. Column density map 
contours are (3, 5, 10, 25, 50) x cr, where the rms noise cr = 1.40 
mJy beam -1 or 2.67 x 10 19 atoms cm -2 . The FWHM synthe- 
sised beam has dimensions 7.6" x 7.6", and is displayed as a 
cross-hatched symbol in the lower left. Middle right: Minor axis 
profile. Bottom left: Major axis profile. 



which was used in the rotation cu rve analysis by lAbe et all 
d!999h and lvan der Kruit et al.l d2001l) . The distance of ESO109- 
G21 is about 30 Mpc distant so the HI beam is relatively large, 
1 179 pc. The HI column density map in Fig.[2blshows the galaxy 
to be well-resolved in the radial direction, while perpendicu- 
lar to the plane the resolution is just sufficient to detect flaring 
in the extreme velocity channel ma ps. Comparison of the XV 
map in Fig.|2b]with that measured by Ivan der Kruit et alj ( 1200 lb 
shows significantly improved structure recovery. However the 
peak signal-to-noise of the HI channel maps is 10.1, which will 
make measurement of the HI flaring quite difficult. 

5.3. ESO115-G021 

This Scd galaxy extends to a maximum radial extent of 445" or 
9.7 kpc (assuming a Hubble flow distance of 4.5 kpc inferred 
from the \'c a i). The scalelength of the projected distribution, hx, 
is 209" or 4.6 kpc, with the edge of the HI disk at 2.lh X - This HI 
column density map (Fig. [3bb displays a very symmetric disk 
which appears to be relatively quiescent with no obvious HI 
extensions above and below the inner disk. The vertical struc- 
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Fig. 5c. ESO146-G014. Top : Peak brightness temperature map. 
Contours are plotted in (60, 80, 100) x cr, where the rms noise 
is 14.9 K. The FWHM synthesised beam has dimensions 7.6" x 
7.6". Middle: Major axis peak brightness temperature profile. 
Bottom: Inferred HI opacity calculated assuming constant HI 
spin temperatures. The resulting maximum opacities along each 
line of sight column through the galaxy disk are plotted for T sp i n 
= 200, 300 and 400 K (bottom to top). 



ture is well-resolved, extending to a height of 1.7 kpc, equiv- 
alent to 8.5 FWHM beamwidths or 2.9h z , where h z is the pro- 
jected scaleheight of 26.7 arcsec or 583 pc . The flaring of the 
vertical HI distribution can clearly be seen in channel maps at 
459,469,479,548,558,568 km s _1 , particularily on the eastern 
side where the flatter domain of the rotation curve is more ex- 
tended (as inferred from the XV map shown in Fig.[3bl. 

The XV diagram does not display barred gas kinematics; 
however a small galaxy like ESOl 15-G021 with a rotation speed 
of 62 km s -1 is unlikely to show the strong orbital patterns of a 
gas bar. The HI column density distribution is dominated by a 
plateau within a radius of 5 kpc, except for a central bright peak 
with a peak brightness of ~ 150 K in the central 1 kpc. 
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Fig. 6a. ESO274-G001. HI column density map overlaid on the 
DSS image. Contours are plotted in (3, 5, 10, 25, 50) x cr, where 
the rms noise cr = 2.02 mJy beam -1 or 2.30 x 10 19 atoms cm -2 . 
The FWHM synthesised beam has dimensions 9.8" x 9.8", and 
is displayed as a cross-hatched symbol in the lower left. 

5.4. ESO138-G014 

ES0138-G14 is a larger Sd galaxy, with a maximum rotation 
speed of 120.4 km s~' and a HI disk extending to 23 kpc. The 
HI column density map in Fig. [4a] shows the HI disk is well- 
resolved by the HI synthesised beam of 10.7" or 960 pc. The 
XV diagram in Fig.l4b1shows depressions in the HI flux at radii 
of ~10 kpc, with the indication of a "figure-8" structure in the 
gas kinematics. This suggests that the galaxy is actually barred, 
as "figure-8" structures are typically seen in gas k inematics of 
barred edge-on systems ([Bureau & Freeman! 1 19991) . 

We were able to observe this galaxy only in the 1.5 A and 
6A ATCA array configurations. Due to the lack of short spac- 
ing observations we are missing much of the extended structure 
on large spatial scales. The channel maps show that the HI disk 
clearly flares vertically with radius, but the missing extended 
structure will hamper accurate measurement of the flaring. 

5.5. ESO146-G014 

This Sd galaxy extends to a maximum radial extent of 118" or 
12.3 kpc (assuming the distance of 21.5 kpc determined from 
its systemic velocity). The galaxy is well-resolved by the syn- 
thesised beam along the major axis, extending to 15 FWHM 
beamwidths on each side of the galactic centre, where each 
FWHM beamwidth is 793 pc. In the vertical direction the mea- 
sured scaleheight of the projected distribution is 2.6 times the 
HWHM beamwidth, indicating that the galaxy is well-resolved 
in the vertical direction also. 

ESO146-G014 has a rotation speed of 70.2 km s" 1 , typical 
of low-mass Sd galaxies. The rotation of ESO146-G014 appears 
solid body from the XV diagram in Fig. [5b] however HI mod- 
elling (see Paper III) shows a differential inner rotation curve 
like other small disk galaxies when observed at high resolution 
(e.g. ISwatersl 119991) . The XV diagram shows no evidence of a 
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Fig. 6b. ESO274-G001. Top left: XV map. XV contours are 
(3,5, 10,25,50, 100) x cr, where the rms noise cr = 2.24 mJy 
beam -1 or 2.56 x 10 19 atoms cirT 2 . The half power beam extent 
over the major axis is shown in the lower left corner. Top right: 
Integrated spectrum. Middle left: HI column density map rotated 
with the galaxy major axis aligned with the X axis. Column den- 
sity map contours are (3,5, 10,25,50, 100) x cr, where the rms 
noise cr = 1.70 mJy beam -1 or 1.93 x 10 19 atoms cm -2 . The 
FWHM synthesised beam has dimensions 9.8" x 9.8", and is 
displayed as a cross-hatched symbol in the lower left. Middle 
right: Minor axis profile. Bottom left: Major axis profile. 



HI bar, which is as ex pected given the low oxygen a bundance of 
6% of the solar value dBergvall & Ronnbacl3.ll 9951) . 

5.6. ESO274-G001 

ESO274-G001 is the closest galaxy in our sample, at a distance 
of 3.4 Mpc in the Centaurus A group. Its proximity allows the HI 
disk to be examined at 160 pc resolution with the 9.8" resolution 
of the HI synthesised beam. To image this galaxy we used two 
partial observations in 1.5 km ATCA array configurations from 
the ATCA archive, and one new 12-hour synthesis observation 
in the ATCA 6A array configuration. From the HI total column 
density map in Fig. [6a] we can see that the vertical axis is very 
well-resolved, although it is clear from the diffuse structure that 
the significant extended HI emission is not observed due to the 
lack of observations in short array configurations. Comparison 
of the HI flux integral with that measured from the HI Parkes 
All-Sky Survey (HIPASS) data shows that about 30% of the HI 
flux is not measured. 

From the XV map in Fig. [6b] we see that the rotation curve 
of this small Sd galaxy rises steeply to a radius of 1 - 2 kpc and 
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Fig. 6c. ESO274-G001. Top : Peak brightness temperature map. 
Contours are plotted in (40, 100) x <x, where the rms noise is 9.4 
K. The FWHM synthesised beam has dimensions 9.8" x 9.8". 
Middle: Major axis peak brightness temperature profile. Bottom: 
Inferred HI opacity calculated assuming constant HI spin tem- 
peratures. The resulting maximum opacities along each line of 
sight column through the galaxy disk are plotted for T sp j„ = 200, 
300 and 400 K (bottom to top). 



then flattens on one side, but continues to rise with a shallower 
slope on the other side. In these outer velocity channel maps, 
where the HI spans a larger range along major axis, it is possible 
to see clear flaring of the gas thickness with radius. However, the 
galaxy is lopsided kinematically, and also slightly lopsided in its 
HI distribution with an indication of a warp on the NE side. 



5.7. ESO435-G025(IC2531) 

ES0435-G25 is the only large spiral galaxy in our sample. It 
is an Sb galaxy with a bright peanut-shaped bulge shown to be 
a bar from the "figure-8" signa t ure in the optical emission line 
kinematics dBureau & Freemanl 1 1999I) . The line splitting faintly 
seen as a figure-8 signature in the HI XV dia gram in Fig.lTblwas 
first detected by iBureau & Freemanl d!997l) . We use the same 
ATCA HI observations in our study. 
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Fig. 7a. ESO435-G025 (IC2531). HI column density map 
overlaid on the DSS image. Contours are plotted in 
(3,5, 10,25,50, 100) x cr, where the rms noise cr = 1.73 mJy 
bearrr'or 2.35xl0 19 atoms cirT 2 . The FWHM synthesised beam 
has dimensions 9.0" x 9.0", and is displayed as a cross-hatched 
symbol in the lower left. 



The HI column density map shows the galaxy extends to 
~25 beamwidths on either side of the galactic centre, and three 
beamwidths on either side of the HI midplane. Above the galac- 
tic centre, the HI column density map in Fig. 1751 shows HI fila- 
ments extending to ~3 kpc above and below the midplane. Like 
ESO109-G021, the other distant galaxy, the peak signal-to-noise 
is quite low at 10.9. Although the HI disk appears to be flared in 
the HI column density map (Fig. l7ab . the flaring is not visibly 
obvious in the HI channel maps due to the low signal-to-noise. 

5.8. UGC7321 

This Sd galaxy extends to a maximum radial extent of 288" 
or 14.0 kpc (assuming the distance of 10 kpc adopted by 
lUson & Matthews! d2003l) l. The scalelength of the projected dis- 
tribution, hx, is 1 12" or 5.4 kpc, with the edge of the HI disk at 
2.6hx- This galaxy exhibits substantial high latitude HI extend- 
ing up to 144" or 7.0 kpc in the inner HI disk. The scaleheight of 
the projected minor axis distribution h z is 1 1 .6 " or 560 pc. Thus 
the high latitude HI clouds extend to a height of 12 scaleheights 
above the equatorial plane. 

The galaxy is well-resolved by the synthesis beam along the 
major axis, extending to 19 FWHM beamwidths on each side of 
the galactic centre, where each FWHM beamwidth is 727 pc. In 
the vertical direction the measured scaleheight of the projected 
distribution is similar to the FWHM beamwidth, suggesting that 
the low latitude vertical structure is dominated by the shape of 
the synthesised beam. However despite the relatively large syn- 
thesised beam, the vertical flaring is clearly visible in most chan- 
nel maps due to the extended roughly flat rotation curve which 
is clearly apparent in the XV map (Fig.[8bl. 

The XV diagram (Fig. [8bl also exhibits the characteris- 
tic "fig ure-8" shape, which can be indicative of barred HI dy- 
namics dKuiiken & MerrifieldL [19951 IBureau & Freemanl 1 19991: 
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Fig. 7b. ESO435-G025 (IC2531). Top left: XV map. XV con- 
tours are (3, 5,10, 20, 30) x <x, where the rms noise cr — 1 .76 mJy 
beam -1 or 2.40 x 10 19 atoms crrT 2 . The half power beam extent 
over the major axis is shown in the lower left corner. Top right: 
Integrated spectrum. Middle left: HI column density map rotated 
with the galaxy major axis aligned with the X axis. Column den- 
sity map contours are (3,5, 10,25,50, 100) x cr, where the rms 
noise cr = 1.68 mJy beam -1 or 2.28 x 10 19 atoms cm -2 . The 
FWHM synthesised beam has dimensions 9.0" x 9.0", and is 
displayed as a cross-hatched symbol in the lower left. Middle 
right: Minor axis profile. Bottom left: Major axis profile. 



lAthanassoula & Bureaulll999HAthanassoulaLl2000l) . Further ev- 
idence for a dynamical bar in UGC7321 was later found by 
Pohlen who detecte d the distinctive boxy-bulge shape in the 
near-IR distribution dPohlen et all 120031) . iPohlen et all (l2003h 
measure the bar length as 1 12 + 21" or 5.4 + 1 .0 kpc. This is in 
agreement with the bar size inferred from the size of the "figure- 
8" kinematics. It also corresponds to the scalelength of the pro- 
jected HI profile. The HI column density has a plateau at radii 
within 5 kpc, declining steeply at radii outside this radius. 

The peak brightness temperature rises steeply at radii inside 
2 kpc to 94 K. Assuming a HI spin temperature of 200 or 400 K, 
the corresponding central HI opacity t is 0.65 or 0.27. The two 
outer peak brightness points correspond to radii of 6 kpc east, 
and 5.5 kpc west, also aligning with the outer Lindblad reso - 
nance points of the stellar bar measured bv lPohlen et all (120031) . 



6. Discussion & Summary 

The resolution of our HI data along the major axis is high, 
with the number of independent beamwidths on each side of 
the galaxy centre ranging from 15 to over 50 in our galaxies. 
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Fig.7c. ESO435-G025 (IC2531). Top : Peak brightness temper- 
ature map. Contours are plotted in (27,46,73) x cr, where the 
rms noise is 9.2 K. The FWHM synthesised beam has dimen- 
sions 9.0" x 9.0". Middle: Major axis peak brightness temper- 
ature profile. Bottom: Inferred HI opacity calculated assuming 
constant HI spin temperatures. The resulting maximum opaci- 
ties along each line of sight column through the galaxy disk are 
plotted for T spin = 200, 300 and 400 K (bottom to top). 



In the vertical direction all galaxies are spanned by at least 5 
beamwidths. Unfortunately three of the galaxies in the sample 
suffer from incomplete imaging, resulting in missing informa- 
tion about the extended spatial structure due to the lack of obser- 
vations along short baselines. For two of these galaxies ESO074- 
G015 and ESO274-G001, the images still contain substantial in- 
formation due to the small linear size of the beam, 290 pc and 
160 pc, respectively. For these galaxies the iterative HI mod- 
elling methods used to measure the deprojected HI density dis- 
tribution and kinematics (Paper II) should provide good mea- 
surements of the flaring. But the HI images of ESO138-G014 
indicate missing extended HI emission and low spatial resolution 
which will probably prevent reliable measurement of the flaring 
in this galaxy. All the other galaxies in our sample are promising 
candidates for accurate measurement of the HI kinematics and 
vertical flaring. 
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Fig. 8a. UGC7321. HI column density map overlaid on the DSS 
image. Contours are plotted in (3,5, 10,25,50, 100,250,500) x 
cr, where the rms noise cr = 0.524 mJy beam -1 or 2.57xl0 18 
atoms cm" 2 . The FWHM synthesised beam has dimensions 
15.0" x 15.0", and is displayed as a cross-hatched symbol in 
the lower left. 



The HI column density distribution of the galaxies in our 
sample varies quite substantially. The four galaxies with max- 
imum rotation speeds i 100 km s" 1 all have HI disks that 
extend to greater than 5 kpc away from the plane. Three of 
them (ES0435-G25, UGC7321 and ES0138-G14) also have the 
"figure-8" signature in the XV diagram sug gestive of an edge-on 
bar in the gas distribution. ESQ435- G25 ([Bureau & FreemanL 
fT999b and UGC7321 dPohlen et al l |2003|) both have boxy- 
peanut shaped stellar bulges consistent with a bar seen edge- 
on. Combined with the figure-8 signatur e in the HI gas struc- 
ture this is strong evidence for a bar dKuiiken & Merrifield , 
1995HBureau & Freeman! 1 19991: 1 Athanassoula & Bureaulll999 ; 



Afhanassoulal 120001) . E138-G014. IC5249 and IC2531 all have 
large HI disks with radii larger than 20 kpc, while UGC7321 is 
much smaller with a radius of only 14.0 kpc. 

The other four galaxies are not much smaller than UGC7321 
in radial extent, but they are considerably smaller in total mass. 
They also lack the vertical extensions in the central disk that sug- 
gest additional sources of heat in the disk. This suggests that one 
of the causes of extended high latitude HI filaments is heating 
related to star formation and other processes associated with the 
bar. 

The peak brightness profiles of each galaxy show that the 
inferred HI opacity (assuming a constant HI spin tempera- 
ture) varies substantially across the major axis of each galaxy. 
The high spatial resolution of our images has made it possi- 
ble to measure high HI brightness temperatures. The maximum 
brightness temperature in our galaxies ranges from 94.4 K for 
UGC7321 to 168.6 K for ES0146-G14. In addition to ES0146- 
G14, two other galaxies have high HI brightness temperatures 
> 150 K. Assuming a mean HI spin temperature of 300 K, the 
maximum inferred opacity for these HI bright galaxies is ~0.7. 
This is comp arable to the maxim um HI opacity of 0.85 measured 
in NGC891 dKregel et all l2004h . Inspection of the peak bright- 
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Fig. 8b. UGC7321. Top left: XV map. XV contours are 
(3,5, 10,20,50, 100,250,500) x cr, where the rms noise cr = 
0.524 mJy beam" 1 or 2.57 x 10 18 atoms cm" 2 . The half 
power beam extent over the major axis is shown in the 
lower left corner. Top right: Integrated spectrum. Middle left: 
HI column density map rotated with the galaxy major axis 
aligned with the X axis. Column density map contours are 
(3,5, 10,20,50, 100,250,500) x cr, where the rms noise cr = 
0.508 mJy beam" 1 or 2.49 x 10 18 atoms cm" 2 . The FWHM syn- 
thesised beam has dimensions 15" x 15", and is displayed as a 
cross-hatched symbol in the lower left. Middle right: Minor axis 
profile. Bottom left: Major axis profile. 

ness as a function of major axis position for these three galax- 
ies shows that, in each case, the region of increased opacity is 
localised spanning a projected radius of ~l-2 kpc. In ESOH5- 
G21, this increased opacity region occurs at the galactic centre. 
However for ES0138-G14 and ES0146-G14, the regions of in- 
creased opacity occur in the outer disk. Three of our eight galax- 
ies appear to be warped, and one is lopsided. This is consistent 
wi th results from lar ger samples, as e.g. mentioned in the review 
bv lSancisi et] (120081) . 
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Fig. 8c. UGC7321. Top : Peak brightness temperature map. 
Contours are plotted in (4, 5,10, 20, 50, 80) x <x, where the rms 
noise is 1.0 K. The FWHM synthesised beam has dimensions 
15" x 15". Middle: Major axis peak brightness temperature pro- 
file. Bottom: Inferred HI opacity calculated assuming constant 
HI spin temperatures of 200, 300 and 400 K. The resulting max- 
imum opacities along each line of sight column through the 
galaxy disk are plotted for T spin = 200, 300, and 400 K (bottom 
to top). 
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Appendix A: Data reduction 

A. 1. Calibration 

Reduction of the HI data was performed with the radio interfer- 
ometry data reduction package miriad. First, the data was loaded 
into miriad with the task atlod, using the barycentric velocity 
reference frame, and specifying the rest frequency of HI. The 
visibilities were then split into separate datasets for each source. 
The primary and secondary calibrators for each source were in- 
spected for the characteristic signs of solar and terrestrial radio 
frequency interference (RFI) using the task uvplt and tvflag. 
Visibilities with RFI, Galaxy HI emission, shadowing or during 
known telescope system problems were flagged bad. The pri- 
mary calibrator data were flux calibrated using the known flux 
with the task mfcal. The secondary calibrator was then used to 
calibrate the telescope gains, and if sufficient parallactic angle 
coverage was obtained with the secondary calibrator it was also 
used to calibrate the bandpass. For observations with low paral- 
lactic angle coverage, such as snapshots, the bandpass solution 
of the primary calibrator was used to calculate the bandpass. 

The calibrated flux solution of both primary and secondary 
calibrators was inspected with uvplt. New secondary calibrators 
were calibrated during observations, and then inspected with the 
task uvflux to see if they resembled an unpolarised point source. 
If the secondary calibrator was too faint, showed polarisation or 
structure in L band then it was replaced by the next nearest suit- 
able bright calibrator for later calibrator scans throughout the ob- 
serving run. Gains and bandpass solutions were inspected with 
the task gpplt. The flux density of the secondary calibrator was 
then corrected using the primary calibrator flux solution. 

The final calibration tables were copied to the target galaxy 
visibility dataset. Before applying the calibration, the gains were 
averaged over the duration of each secondary calibrator pointing 
(3 minutes) and gains interpolation was limited to the interval 
between secondary calibrator observations (30-45 minutes). 

The calibrated target galaxy visibilities were then inspected 
with the tasks tvflag, uvplt, and any data contaminated by RFI 
or shadowing was flagged bad with the appropriate flagging tool 
(tvflag, uvflag or blflag). 

A.2. Continuum subtraction 

To identify the presence of radio continuum sources in the field 
of view, a large 30' channel-averaged (chO) image was made us- 
ing the task invert extending out to approximately the half power 
radius of the primary beam. The positions and flux in radius were 
fitted for any continuum sources found in the chO map using the 
task cgcurs. The visibility spectra of each baseline were then 
inspected using the task uvspec Spectral channels free of HI 
emission were identified. Observations of nearby galaxies often 
showed Galactic HI emission within the observed bandpass. It 
was not necessary to flag these channels. Instead they were ex- 
cluded from the list of line-free channels used in the continuum 
subtraction. As the target galaxy HI emission was always sepa- 
rated from the Galaxy line emission by at least 1 00 - 200 km s~ 1 , 
they were also not included in the galaxy HI image cubes. 

The radio continuum was then removed by performing a low 
order polynomial fit to the line-free channels with the task uvlin. 
If the observation included a radio continuum source with flux 
significantly brighter than the channel-averaged HI flux of the 
target galaxy, then the phase center was shifted to the position 
of the brightest contaminating continuum source during the con- 
tinuum subtraction. The continuum-corrected visibility spectra 



were checked with the task uvspec and a new large chO image 
was formed to search for additional fainter residual continuum 
sources that remained after the first continuum subtraction. If 
found, continuum subtraction was repeated with the phase cen- 
ter shifted to the next brightest source. The process was repeated 
until the chO image showed no significant contaminating contin- 
uum sources. For observations that were taken during daylight 
hours, the continuum subtraction was performed with the phase 
center shifted to the position of the Sun. This procedure could 
be compromised by continuum emission from the galaxy itself 
or background sources and in all cases it worked satisfactorily. 

For the narrow bandpass VLA observations, it was not pos- 
sible to span the full galaxy emission in a single bandpass when 
using the high spectral resolution correlator configuration. In 
such cases, the galaxy was imaged with two overlapping band- 
passes. As these bandpasses included only a few line-free chan- 
nels on one side of the bandpass it was not possible to perform an 
accurate 1st order polynomial fit to the line-free channels with 
uvlin. Instead any bright continuum sources were subtracted 
from each bandpass with the task uvmodel using a point source 
model with the fitted flux and position of the radio source. Any 
residual continuum emission was removed with a zero-order fit 
to the few line-free channels in each bandpass. This performed 
a satisfactory continuum subtraction for the three galaxies ob- 
served at the VLA. 



Appendix B: Imaging 

A preliminary dirty HI data cube and dirty beam map was 
formed using the default gridding where each pixel size is 1/3 
of the synthesised beam FWHM along the X and Y axes of the 
beam. The accurate beam shape was then measured from the 
central peak of the dirty beam by fitting an elliptical Gaussian to 
determine the major 9„ W j and minor mi „ axes and position angle 
of the synthesised beam. The preliminary data cube was used to 
determine the size and velocity range of the HI emission. A new 
dirty datacube was then made with square pixels each of 1 /3 or 
1 /5 of the synthesised beam FWHM major axis 9 ma j. 

The visibilities were weighted using a robust weighting of 
0.4 dBriggsL 1 19951) to provide a optimal compromise between 
spatial resolution and suppression of sidelobes. A positive robust 
weighting causes the visibilities in the more sparsely-sampled 
parts of the visibility phase space to be down-weighted relative 
to uniform weighting. This reduces the amplitude of the side- 
lobes, without significantly increasing the size of the synthesised 
beam FWHM. The barycentric velocity reference frame and the 
radio definition of velocity was used for the velocity axis of the 
HI data cube. We chose the radio definition because it yields a 
constant velocity increment between channels, which was useful 
for later modelling of the HI distributionQ 



9 There are two common approximations to the relativistic equation 
describing the doppler effect - the "radio definition" for which equal 
increments in frequency translate to equal increments in velocity, 

v 

Vradio = C(l ), 

vo 

and the "optical definition", 

^optical — C\ I)- 

V 

The difference in v for fixed v is ~ 1 % for 2500 km s~ 1 , the highest radial 
velocity in our sample. 
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We used the miriad clean algorithm to correct for effects in- P.C. van der Kruit or K.C. Freeman at the email address given 
traduced by the sparcity of the visibility sampling. To find the on the first page, 
approximate HI emission region in a each channel, first an ini- 
tial clean was performed over the central quarter of each chan- 
nel map, and a low resolution data cube was made by restoring 
the clean components with a Gaussian of twice the full resolu- 
tion synthesised beam size. The noise in each plane of the low 
resolution cube was then measured, and inspected to find the re- 
gion of HI emission in each plane. Following standard practise, 
a clean mask was created by selecting the region in each plane 
of the low resolution data cube with HI emission greater than 
twice the rms noise (2<x) using the task cgcurs. The noise level 
of the full resolution dirty cube was measured and then the ar- 
eas defined by the clean mask were cleaned down to 0.5cr using 
miriad clean. We realise that this is not the optimal procedure 
to find the low level emission as CLEAN components found in 
a high resolution image usually do not represent the low level 
extended emission properly. Better practice would have been to 
make low resolution images from the continuum subtracted u,v 
data directly and CLEAN those. In the case of this study, how- 
ever, only high resolution images are produced and used. 

The final full resolution HI data cube was produced using 
the task restor which produces a residual map for each channel 
map by subtracting the transformed clean components, and then 
adds the restored flux to each residual image by convolving the 
clean components by the clean beam. A circular Gaussian was 
used for the clean beam, with the FWHM set by the measured 
major axis FWHM of the dirty beam. After restoration the full 
cube was corrected for the effect of the antennae primary beam 
using the task linmos. The size of the restoring beam is shown in 
column 2 & 3 of Table [2] 

As multiple observations in different array configurations 
were combined to provide the maximum visibility sampling, res- 
olution and sensitivity, some galaxies were observed at slightly 
different pointing positions. These observations were imaged us- 
ing the mosaic mode in invert, and deconvolved using the task 
mosmem which performs a maximum entropy deconvolution of 
data with multiple pointings. 

In the case of VLA data for which multiple bandpasses were 
needed to span the galaxy HI emission, the visibility data were 
split into separate bandpasses prior to imaging, with the over- 
lap region combined using the task uvcat. Each bandpass was 
then imaged, cleaned and restored separately, before concatenat- 
ing the final subcubes with the task imcat. Although the channel 
separation at the joins was not exactly the same as the channel 
width elsewhere, the difference was less than 1 km s _1 , which 
is unlikely to effect subsequent modelling of the HI distribution. 
The overlap channels were observed for a longer integration time 
and consequently these channel maps have a lower rms noise. 
However as channel maps and moment analysis were made us- 
ing the noise measured in each channel map, the difference in 
noise level over the velocity range of the galaxy is handled ap- 
propriately and does not cause noise contamination in the de- 
rived images. The measured noise levels of the HI data cubes are 
shown in Table [3] 

After forming the full resolution HI data cubes with Miriad, 
all subsequent imaging, moment analysis and measurements 
were undertaken using the Interactive Data Language (IDL). 

Appendix C: Channel maps 

In (lQ'BrienLl2007h all the channel maps that contain HI emission 
have been dispayed as contour diagrams. Interested readers can 
obtain an electronic .pdf file with these maps upon request from 



